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Abstract 

The momentum-dependence of the p^-uj mixing contribution to charge- 
symmetry breaking (CSB) in the nucleon-nucleon interaction is compared in 
a variety of models. We focus in particular on the role that the structure 
of the quark propagator plays in the predicted behaviour of the p^-uo mixing 
amplitude. We present new results for a confining (entire) quark propagator 
and for typical propagators arising from explicit numerical solutions of quark 
Dyson- Schwinger equations. We compare these to hadronic and free quark cal- 
culations. The implications for our current understanding of CSB experiments 
is discussed. 
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A wide variety of charge symmetry breaking (CSB) phenomena in nuclear 
physics appear to be well explained in terms of one-boson exchange poten- 
tials, when electromagnetic effects, the neutron-proton mass difference, and 
isoscalar-isovector meson mixing are included. A recent and extensive review 
of CSB has recently appeared||l|. Ultimately we would like to understand CSB 
in terms of electromagnetic effects and the u-d quark mass difference in a mi- 
croscopic quantum chromo dynamics (QCD) based description of the strong 
interactions. Any potential which gives rise to CSB in the nucleon-nucleon 
(NN) sytem can be categorized as either a class III or class IV interaction in 
the conventions of Henley and Miller 0. 

In the relatively successful explanation of a number of observed CSB phe- 
nomena, calculated contributions from p^-uj mixing play a major role. These 
phenomena include: (i) differences between pp and nn scattering lengths and 
effective range parameters, (ii) binding energy difference of and ^He, (iii) 
analyzing power difference in n-p elastic scattering, and (iv) the Nolen-Schiffer 
anomaly. In the one-boson exchange model of the nucleon-nucleon interaction 
the p^-uj potential is generated by p^-u mixing in the intermediate vector- 
meson propagators!^. Another source of CSB arising from meson mixing is 
that from vr^-r] mixing, but this is typically much less important than that due 
to the p° and uj. The long standing problem of explaining the binding energy 
differences of mirror nuclei, (the Nolen-Schiffer anomaly) 0], appears to have 
been largely resolved as a result of the p^-u mixing potential 0. The theoret- 
ical understanding of the difficult neutron-proton analyzing power difference 
measurements carried out at TRIUMF0 and IUCF||^ relied on p^-u mixing 
and in particular this was the dominant contribution in the latter. [|], ^ p!0| . 

Although the situation appears to be very satisfactory, there is however 
at least one potential problem with the standard evaluations which has been 
pointed out recently by Goldman, Henderson and Thomas(GHT)|]rT|. The 
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problem is associated with the assumed off-shell behavior of the p^-uo mixing 
matrix element. The standard assumption is that the mixing amplitude is mo- 
mentum independent with its value extracted from the experimental data on 
e^e~ TT+TT^ at the uj pole, ^ m^, while the exchanged mesons in Eq. @ 
have space-like four-momentum, < 0, and are therefore highly virtual. In 
their initial study a simple model was used where the vector mesons are con- 
sidered as quark-antiquark composites. Because of the the up-down quark 
mass difference the p^-uj mixing amplitude is generated by an intermediate 
quark loop. A p° of four momentum dissociates into a quark-antiquark pair 
and then recombines to form an cu, or vice versa. The quark and antiquark 
momenta are C.^ + and —i^ + {q^ 1'^) respectively, where C. is the loop 

momentum. There will also be a contribution from electromagnetic effects, but 
this will be much smaller than that from the quark mass difference. Simple ver- 
tex functions were introduced and free quark propagators were used. Once the 
mixing amplitude was calculated, GHT obtained the coordinate space (static) 
potential by Fourier transforming the resulting momentum-space potential. 
The basic conclusions of the GHT calculation were: (i) the mixing amplitude 
is strongly momentum dependent, (ii) there is a node in the potential at about 
0.9fm and, (iii) because of the node, the potential changes sign and so its 
importance is greatly suppressed. 

The potential implications of this result of GHT are profound in view of 
the central role meson mixing plays in our understanding of charge symmetry 
breaking. If these results can be confirmed in other more realistic calculations, 
it may be necessary to find new sources of CSB in order to explain existing data. 
Since the initial investigation by GHT other studies appear to support their 
basic conclusions. These include investigations of tt^-// mixing using a quark 
loop model[|T^, chiral perturbation theory[^], and a hadronic model|n|. In 
addition, each of two further studies of p^-uj mixing using a hadronic model ]TB[] 



3 



and using QCD sum rules|T6| also support the basic conclusions of GHT. 

In the present study we investigate the role of quark confinement and the 
nature of the quark propagator in a q-q based description of the p^-uo mixing 
amplitude. This study is motivated by one of the difficulties associated with 
the GHT calculation, namely that, as their work did not include quark con- 
finement, an unphysical gg-pair production threshold resulted in the timelike 
region - at = 4mg. This corresponds to the vector meson being able to 
dissociate into an on-shell quark-antiquark pair for q^ > 4m^. We review some 
necessary formalism before discussing details of our quark propagators. 

The interaction Lagrangian densities for the nucleon-p and nucleon-u; cou- 
plings are 

d 

Cnnp = QpipY^ ■ i^Pt^ + fpi^(^'"'r ■ V'^Pv , (1) 
'Cnnoj = Quj^Yi'^p ■ (2) 

These definitions are standard, (although some authors vary in their definitions 
of the couplings by including additional factors of 1/2). No tensor coupling 
for the NNuo vertex is included since we choose typical coupling constants 
and form factors determined by the Bonn group |T^. The contribution from 
p^-uj mixing to the A^A^ potential, is given by (with = i(j^"qy/2Misi and 
Cp = fp/gp) 

V.TM = ng)7^(l)7.(2)[r.(l)+r.(2)] , (3) 
V,7{<!) = l^(g)C,[r^(l)7,(2K(l)-7^(l)r,(2)r.(2)] , (4) 

where 

X _ 9p9c.FNpiq^)FNUq^){p°\H\ uj) 

where {uj\H\p'^) is the meson mixing amplitude as defined in Ref. [Q]. We 
have introduced the standard Bonn form factors at the nucleon-meson vertices, 
which for on-shell nucleons are given by FArp(g^) = — (g^/A^)] andF^(g^) = 
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1/[1 — (g^/A^)]. With this normahzation for the form factors, the appropriate 
Bonn couphngs are gl/A-n = 0.41 and gl/4:TT = 10.6. Also, Ap = 1400 MeV, 

= 1500 MeV, Cp = fp/gp = 6.1 and as aheady stated C^^ = fui/gui = 0. 

Our starting point for the present calculations is the expression for the puj 
mixing self-energy in terms of a quark-antiquark loop insertion 



4/ 



u — > d 



n^'^(g^) = ?>ig,pg,.j ^ tr[r^(£+,£_)5„(/_)C(£_,£+)5„(£^ 

(6) 

The trace is over spinor indices and the factor of three comes having per- 
formed the colour trace. T^{p\p) and r^(p',p) are the quark-vector-meson 
vertex functions for the p° and u respectively, Si{q) is the quark propagator 
with flavour i {i = u,d) and i± = i ± (q'/2). We take for simplicity here 
^p{p',p) = l''Fqp{p',p) and T^{p',p) = 'y^'Fg^{p' ,p), where Fgp and Fg^ are the 
quark-vector-meson form factors. Since we wish to study the g^-dependence 



arising from the quark propagator in the loop, we follow Ref. [11| and con- 
sider these form factors to be functions of £^ only, i.e., Fgp(£+,£_) = Fgp^i"^) 
and Fgo; = We do not wish to introduce any unconstrained 

g^-dependent behaviour through the quark-meson vertex functions. 

We can make use of Lorentz invariance and current conservation (i.e., 
q^U^"' = 0) to write W = [ - g^"" + {q^q'" / q^)]Ii{q^), where this is the def- 
inition of the scalar function n(g^). Note that in the case of the photon 
self-energy we would require n(0) = in order that the photon not acquire 
a mass through self-energy insertions, however there is no such restriction in 
the present case. Since we are only interested in coupling to conserved exter- 
nal currents, (i.e., since q^J^ = for the nucleons), we need only retain the 
g'^'^Il{q'^) part of Eq. (P) in our calculation of the A^A^-interaction. A com- 
parison of the definition [Q] of the mixing amplitude {uj\H\p^) with the usual 
Feyman rules shows that this is identical with the pu mixing self-energy n(g^), 
i.e., {uj\H\p^) = n(g2). 
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The general form of the quark propagator in Minkowski space is 



1 



(7) 



A{q^)4-B{q^) q^-M^{q^) 
where A{q'^) and B{q^) are scalar functions of q^ and, Z{q^) = 1/A{q'^) and 
M{q^) = B{q^)/ A{q^). One possible mechanism of quark confinement is 
that the quark propagator does not have a mass pole, i.e., that the function 
= Z{q^)/[q^-M^{q^)] does not have a pole or is an entire function in the 
complex q^ plane [|T8[] . One sees that there is no mass pole when the function 
Z{q^) goes to zero for q^ — > M'^{q^). One explicit quark model where this 
occurs can be found in Ref ([l^) where this property arises from the solution 
of a model quark Dyson-Schwinger equation. 

We can isolate the g^'^'n.{q^) part of 11^^ and we then find in our gg-model 
that the mixing amplitude can be written as 



{qJl 
q2 



Ae + q^ + 4:Mu{f)Mu{il 



u 



d 



where for the w-quark propagator Su{p) = -Fu(p^)[^+^«(p^)] and similarly for 
the (i-quark. The quark- vector-meson couplings are estimated from a standard 
quark model analysis of the nucleon-vector-meson coupling. We find that 
Qqp ~ (3mg/5MAr)5(p(l + Cp) and g^^ ~ {mq/MN)gp{l + C^) c:^ {mq/MN)gp. 
With quark masses of ~ 400MeV this gives g^p ^ 3.9 — 4.3 and gq^^ ~ 4.8 — 5.2 
as typical ranges for these couplings. Throughout this work we will use then 
gqpgqi, = 20.0. 

In the chiral limit, (i.e., with massless quarks and neglecting electromag- 
netic and weak effects), the pion decay constant can be related to the quark 



propagator through the Goldberger-Treiman relation (see, e.g., [20| and refer- 
ences therein) to give 



f 



A2i) 



dH F{f)M{f 



(27r)4 £2 _ M2(£2) 



2^ de^ ^ ' 



(9) 
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Factors of three and four follow from colour and spinor traces respectively. 
(Note that this differs slightly from the expression in Ref. which used 
a dressed axial-vector vertex instead of the appropriate bare one. The bare 
vertex avoids double counting problems). Similarly, in the chiral limit (see, e.g., 
pOU and references therein) we obtain an expression for the quark condensate 
in terms of the quark propagator 

(qq) = = Tr5(£) = (-12.) J 0^,F{f)M{f) , (10) 

where the trace is over colour and spinor indices. We use these two results to 
constrain the quark propagators by fitting them to the experimental values of 
U = 93MeV and {qq) = = -(225MeV)^ (i.e, mo = 225MeV). 

Motivated by the confining quark model from the Dyson- Schwinger equa- 
tion study of Ref. |T^ and the need to avoid the artificial quark pair-production 
threshold encountered by GHT in the timelike region we first consider the fol- 
lowing form for the quark propagator [c.f., Eq. (|^)]: 

F(g2) = foe-'^'/^'' and M{q^) = m , (11) 

where /o and /x are parameters and m is to be chosen as a typical infrared quark 
mass. This propagator is confining and analytic everywhere in the complex 
q^ plane (i.e., entire). It captures the essential elements of the propagator 
of Ref. [0 and has the great advantage that with the appropriate choice of 
gaussian quark-vector-meson form factors we can use Eq. (P) to obtain an 
analytic expression for the mixing amplitude n(g^) in both the timelike and 
spacelike regimes. In practice it is easier to calculate the expressions for /^^ 
and the condensate mo [i.e., Eqs. (§,0)] after first rotating to Euclidean space. 
Evaluating these expressions for our analytic (confining) propagator gives 



/o = 47r2mi^/3m/ 

(12) 



2 ^ 3/omV / „ /-^ .. e t 

J -IT 



47r2 



ae" I dt ^ 



a={rn/ fi)'^ 



Hence we can now fix the parameters /o and jj,. Motivated by the infrared 
values of the dynamically generated quark mass in Dyson- Schwinger equation 
studies [^] we choose here m = 450MeV for the "constituent" quark mass. 
Then substituting the expression for /q into the expression for in Eqs. ( [12D 
lets us solve for /i to give ^ = 523.6MeV and hence we then find /o = 4.432 x 
10~^MeV^^. Thus all parameters are now constrained. 

Since the u-d quark mass difference is commonly estimated to be ~ 4MeV 
we use md—rriu = 4MeV throughout this work. Hence the "constituent" masses 
are m„ = 450MeV and = 454MeV and fou = fo and fi are as determined 
above. The simplest and most reasonable assumption is that fi has the same 
value for the u and d quarks, and hence from Eq. (|12]) fou/ fod = rnd/niu. 
As discussed above in order to facilitate the integration in Eq. (H), we use 
Fqpi.^"^) = Fqi^i.^^) = exp(£^/A^) with A = IGeV, which is a typically reasonable 
number for the form factor cut-off. We see that these form factors fall off as 
gaussians in the spacelike region. The integral in Eq. can be performed 
easily and one obtains for H(g^) 



/oJ V - 4m^ - (4/A)] 



u d 



(13) 



167r2A2 

where A = 2[(1/A2) + (l/zi^)]. The mixing amplitude has thus been calculated 
for all timelike and spacelike for the analytic (confining) case. 

We can also numerically evaluate H(g2) for numerical solutions of quark 
Dyson-Schwinger equations. These solutions are fitted to both f^^ and the 
quark condensate in the usual way||20|| and in addition have the correct per- 
turbative asymptotic behaviour in the spacelike (Euclidean) region. Unfortu- 
nately it is not yet known how to continue these solutions into the timelike 
regime in the general case. We have introduced a 4MeV mass splitting of 
the infrared mass and have solved for the u and d propagators in the usual 



waypm. The quark- vector-meson couplings {gqp and g^^) have already been 
deduced and for the corresponding form-factors we choose the form used by 
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GHT for ease of comparison, i.e., Fqp{e) = F^^{e^) = 1/[1 - (f/A^)]. We 
again use A = IGeV. We can now use Eq. to evaluate n(g^) for a typical 
Dyson-Schwinger equation solution in the spacelike region. 

Finally, for comparison we repeat the calculation with the numerical Dyson- 
Schwinger solution replaced by a free quark propagator with artificially heavy 
quarks (m^ = 600MeV and = 604MeV). This calculation is then identical 
to that of GHT except that now there is no troublesome quark threshold 
between the spacelike region and the a;-mass point where the mixing amplitude 
is actually measured. 

In order to understand the importance of the momentum dependence of 
the mixing amplitude on the nucleon-nucleon potential due to vector meson 
exchange, we consider the static, central part of this potential. Besides a 
central potential, the exchange of a vector boson gives rise to spin-orbit, spin- 
spin and tensor potentials. Although the spin- dependent potentials are the 
relevant ones for the TRIUMF and lUCF experiments, the radial part of these 
potentials will be directly related to that of the central one. The coordinate 
space A^A^ potential due to pcu-mixing, V{r), is the Fourier transform of V{q) 
in Eq. @, where g° = and where we here define q = \q\ 



In Fig. 1 we show the momentum-dependence of the p^-uj mixing ampli- 
tude for the three quark calculations described here, i.e., the confining (an- 
alytic), Dyson-Schwinger equation (dse), and heavy- free (free) cases. Also 



tially parameter free calculation was made using an A^A^ loop and the n-p 
mass difference. The two hadron-model curves differ in that one (hadronJf) 
has the usual Bonn form factor applied at the vertices of the nucleon loop in 




(14) 



shown for comparison are recent hadronic calculations |]T5[, where an essen- 
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the spaceklike region as well as at the meson-nucleon vertices. We see that 
all calculations appear capable of fitting the data with typical parameters and 
that the momentum dependence is remarkably similar in each case. We see 
strong momentum-dependence of the mixing amplitude and a change in sign 
in the vicinity of = 0. This is also what was seen in the recent QCD sum 
rules study although there Il{q^) was seen to change sign in the timelike 
region at ^ 0.25 — 0.40GeV^ rather than the spacelike region. 

In Fig. 2 we show the momentum-space potential V{q) given in Eq. (|^) for 
the various models. Also for comparison we show the result obtained when 
the usual assumption of momentum-independence of the mixing amplitude is 
made. In Fig. 3a) we show the corresponding coordinate-space potentials V{r) 
given by Eq. (0) for each of these cases. Fig. 3b) has an enhanced vertical scale 
in order to show the location of nodes in the potential V{r). We see that all 
calculations predict strong suppression of the p^-u mixing potential compared 
with that resulting from the usual momentum-independent assumption for 
the mixing amplitude. We see that the Dyson- Schwinger propagator and the 
hadronic models predict the opposite sign at small r and have nodes between 
0.5 and 0.8fm, whereas the heavy free quark and confining (analytic) quark 
cases have the same sign as the usual assumed potential, but are very strongly 
suppressed. It is interesting to note that the heavy free quark case presented 
here is identical to the case studied by GHT with the exception of having 
heavier (i.e., 600MeV) quarks. We see that this increase in mass has removed 
the node in V{r) for the free propagator case. Given the variety of approaches 
to the calculations, the results are remarkably similar. 

The possible implications for our understanding of CSB are far-reaching. 
For example, in the case of the p-u mixing contribution to the class IV CSB in 
n-p elastic scattering ||^-[rU| there is a competition between large-r fall-off in 
the potential and short- distance suppression of the distorted A^-A^ wavefunc- 
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tion. The result is that the conventional contribution peaks around 0.9fm. It 
is interesting that this is the region where nodes occur. All models predict 
at least strong suppression in this region. It is clear that the previously as- 
sumed theoretical understanding of charge-symmetry breaking phenomena is 
now questionable. The models and treatments to date appear to be in over- 
all agreement. The crucial question is whether these conclusions will survive 
future, more elaborate examinations. 
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Fig. 1 The p-uj mixing amplitudes calculated in a variety of models, (see 
text). Also shown for comparison is the experimentally extracted amplitude 
at the a;-mass-shell point and the usual assumed g^-independent behaviour. 

Fig. 2 The p-uj momentum space potentials, V{q), for the mixing amplitudes 
shown in Fig. 1. V{q) is defined in the text. 

Fig. 3 The p-uj coordinate space potentials, V{r)^ for the mixing amplitudes 
shown in Fig. 1. V{r) is defined in the text. We show the full range of the 
potentials in a), while in b) an enhanced vertical scale is used to show in detail 
how the potential changes sign in some cases. 
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